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Long-lived population storage in optically pumped levels of rare-earth ions doped into solids, re¬ 
ferred to as persistent spectral hole burning, is of significant fundamental and technological interest. 
However, the demonstration of deep and persistent holes in rare-earth ion doped amorphous hosts, 
e.g. glasses, has remained an open challenge since many decades - a fact that motivates our work to¬ 
wards a better understanding of the interaction between impurities and vibrational modes in glasses. 
Here we report the first observation and detailed characterization of such holes in an erbium-doped 
silica glass fiber cooled to below 1 K. We demonstrate population storage in electronic Zeeman- 
sublevels of the erbium ground state with lifetimes up to 30 seconds and 80% spin polarization. In 
addition to its fundamental aspect, our investigation reveals a potential technological application of 
rare-earth ion doped amorphous materials, including at telecommunication wavelength. 


Persistent spectral hole burning in cryogenically cooled 
rare-earth ion (REI) doped solids, the process of trans¬ 
ferring atomic population via optical pumping from the 
ground level of the REI into a long-lived (metastable) 
level, has been extensively studied for five decades mm- 
The observation of persistent holes is of fundamental in¬ 
terest as they allow characterizing the interaction be¬ 
tween REIs and their environment, and they find appli¬ 
cations in both classical and quantum information pro¬ 
cessing [3]. Deep, narrow, and long-lived holes have been 
observed for various REIs doped into crystalline hosts 
HHS]. These holes generally result from population stor¬ 
age in nuclear spin hyperfine levels, or, for REIs having 
an unpaired electron spin (so-called Kramers ions), in 
electronic Zeeman sublevels. Compared to REI-doped 
crystals, REIs in glasses experience much larger inhomo¬ 
geneity of optical and spin transitions, as well as coupling 
with inelastic tunneling modes (also referred-to as two- 
level systems - TLS) [THS] that results in significantly dif¬ 
ferent decoherence and population relaxation properties. 
Persistent population transfer due to coupling with TLS, 
photochemical hole burning, and using spin states have 
been reported more than two decades ago [TOUI^ . but 
the interaction between REIs and vibrational modes in 
amorphous media - the generalization of the well-known 
phonons that are present in crystals - is still not well 
understood. Furthermore, the observed holes have been 
shallow or hundreds of MHz wide, making them unlikely 
to be suitable for information processing applications. 

An exception is the erbium-doped fiber that we have 
used recently for quantum state storage [H HI- The 
goal of the present paper is to detail its hole-burning 
mechanism. We observe deep, persistent holes with 15 
MHz width (limited by power-broadening) and lifetimes 
approaching a minute that arise from population redis¬ 
tribution among the Zeeman-split electronic sublevels of 
the Er^"*" ground state. This improves upon the prop¬ 


erties observed before for REI-doped glasses in terms of 
spin polarization (i.e. hole depth) and spectral width 
[H [H]. We furthermore characterize the dependence 
of the spin relaxation rate on magnetic held, tempera¬ 
ture, wavelength, and erbium ion concentration. This 
allows us to determine the processes responsible for elec¬ 
tron spin state relaxation in this amorphous medium and 
provides the necessary information to optimize operation 
parameters for practical applications. 

Except where mentioned otherwise, our experiments 
employ a commercially available, 20 m-long silica hber 
co-doped with 190 ppm erbium (INO S/N 404-28252). 
The hber is cooled to temperatures T as low as 0.65 K 
using an adiabatic demagnetization refrigerator, and a 
magnetic held B of up to 0.25 T that lifts the Kramers 
degeneracy of the ground and excited electronic levels, 
splitting each into two Zeeman sublevels. We burn spec¬ 
tral holes into the erbium ions’ inhomogeneously broad¬ 
ened ‘^/i 5/2 Iiz /2 zero phonon line using 500 ms long 
laser pulses with a peak power of 1-5 /rW derived from a 
narrowband, intensity-modulated, continuous wave laser 
operating, except where otherwise specihed, at 1532 nm 
wavelength. After a variable waiting time that exceeds 
the 11 ms excited state lifetime, we probe the absorption 
prohle by linearly chirped the frequency of the laser light 
over 500 MHz across the generated hole. The 1 ms-long 
scans are implemented using a phase modulator driven 
by a serrodyne signal |16j . Depending on the experimen¬ 
tal conditions, we observe persistent spectral holes with 
line widths as narrow as 15 MHz, including the effect 
of power broadening. However, from coherence lifetime 
measurements, we know that homogeneous linewidths as 
narrow as a few MHz can be achieved in our fiber m- By 
measuring the area of the spectral hole with respect to 
the waiting time, we characterize the decay dynamics of 
the spectral hole and extract the underlying population 
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FIG. 1: Spectral hole area as a function of the waiting time for 
B = 500 G and T — 0.65 K. The experimental data is fitted by 
the sum of two exponential decays with characteristic times 
Ta = 1.32 ± 0.09 s and Th = 25.9 ± 1.5 s. The inset shows the 
fiber absorption profile at T = 0.8 K (dashed blue line) and 
the hole lifetime Tf, for various wavelengths (black dots). 

relaxation processes. 

First, we determine the decay dynamics for persistent 
holes, with a typical example shown in Fig. We find 
that all decays, for temperatures ranging from 0.65 to 3.5 
K and magnetic fields from 0 to 0.25 T, can be described 
by two exponential functions with similar weights and de¬ 
cay times {To) on the order of a second and (T),) on the 
order of a few tens of seconds. This observation could in¬ 
dicate the presence of two different classes of erbium ions 
in the fiber, or could result from the average over a non- 
uniform, continuous distribution of decays corresponding 
to ions with different spin-state lifetimes. In addition, we 
characterize the wavelength dependence of the persistent 
hole lifetimes at a temperature of 0.8 K; see inset of 
Fig.a We observe long-lived spectral holes with compa¬ 
rable lifetimes across the entire THz-wide inhomogenous 
absorption line. 

Second, we assess the efficiency of the hole burning 
mechanism, which we define as the maximum degree of 
spin polarization. It characterizes how deep a persistent 
spectral hole can be burnt. Spectral holes with absorp¬ 
tion reduced by as much as 80% compared to the opti¬ 
cal depth of the erbium-doped hber before spectral hole 
burning were observed at T = 0.7 K, and with a probing 
delay of 50 ms. Furthermore, as the temperature rises 
from 0.7 to 3.5 K, the hole depth decreases linearly. In 
consequence, we extrapolate that it is, at least in prin¬ 
ciple, possible to burn a 100% deep hole at an optimal 
temperature of 0 K. This fact alone excludes population 
redistribution via the TLS mechanism, as the latter gen¬ 
erally only allows for persistent holes with a maximum 
depth of 50%, and typically much less - even at T = 0 K 
m- Furthermore, given the width of the observed holes 
and the small intensity used to create them, we can also 
exclude photochemical hole burning. 

An unambiguous way to conhrm our hypothesis of Zee¬ 



FIG. 2: Magnetic field dependent depth of a 200 MHz-wide 
hole normalized to the depths of a 100 MHz-wide hole at each 
field, which is not significantly affected by population transfer 
from an anti-hole. Please note that the absolute hole depth 
of a narrow hole decreases at B > 0.06 T due to decreasing 
lifetimes, as shown in Fig. [^(c). Experimental data is shown 
by black dots; simulated results are in red. 

man level storage as the hole burning mechanism would 
be to burn one narrow spectral hole, observe the re¬ 
gions of increased absorption due to population redistri¬ 
bution, i.e. anti-holes, and extract their magnetic field- 
dependent shift. However, this was not possible in our 
case due to the disorder in the amorphous material, which 
leads to very broad anti-holes that cannot be resolved. 
To indirectly verify the presence and positions of these 
anti-holes, we employ a novel method that is based on 
burning a wide hole and measuring its depth as a func¬ 
tion of magnetic field. The results for a 200 MHz-wide 
hole are shown in Fig. As the field decreases, the 
anti-holes begin to overlap with the central hole, causing 
it to become more shallow (due to the first-order elec¬ 
tronic Zeeman splitting, the frequency shift between the 
hole and anti-hole depends linearly on the magnetic field 
strength). Modeling the hole and the anti-holes with 
Lorentzians of magnetic field-dependent width and po¬ 
sition, we can predict the change in hole-depth as a func¬ 
tion of the magnetic field, see Fig.[^ This allows estimat¬ 
ing a mean hole/anti-hole splitting of 25 GHz/T, which 
is comparable to the values typically measured in erbium 
doped crystals [MUD]. This further supports our con¬ 
clusion of population storage in Zeeman levels. We also 
find that the anti-hole width is about six times larger 
than the splitting between hole and anti-hole, which re¬ 
flects the large inhomogenous broadening of the Zeeman 
transition in the fiber, as discussed earlier. 

To gain insight into the relaxation mechanisms that 
determine the observed spin lifetimes as well as to iden¬ 
tify optimal conditions for potential applications based 
on persistent spectral hole burning, we experimentally 
and theoretically study magnetic field and temperature 
dependence of hole decay rates (the inverse of the hole 
lifetime); examples of the experimental data are shown in 
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Figs. [^and|^ We note that the two decay components 
Ta and Tf, exhibit very similar behavior as a function of 
magnetic field and temperature for the assessed range of 
parameters. Given the difficulty of measuring the longer 
Zeeman lifetime with high precision at high temper¬ 
ature, we consider only the fast decay component Ta in 
the following analysis. 

To describe our experimental data, we develop a 
model based on the framework proposed in m for an 
Er:Y 2 Si 05 crystal, but adapt it to an amorphous host, in 
which the Debye model of phonons is not applicable and 
has to be replaced by TLS and local vibrational modes. 
We assume the spin relaxation rate to be described by 
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TABLE I: Parameters for Eq. [^to fit the experimental data, 
with Oi in 10® s“®, 0:2 in Hz K~^ and 03 in Hz 


1 

Ta/b 


ai 

T% + ^B 


sech^ 


/ (m}B\ 
2kT ) 


+ a2B^T'^ + ■ 


( 1 ) 

The first term corresponds to the average mutual spin 
flip-flop rate due to magnetic dipole-dipole interactions 
between erbium ions, with g the (/-factor of the Er ions, 
Hb the Bohr magneton and k the Boltzmann constant. 
We assume g = 9, which is consistent with averaging over 
all possible direction-dependent values, ranging from 0 to 
18 P^27) . Provided the erbium ions are uniformly dis¬ 
tributed, the coefficient ai scales quadratically with Er 
concentration (we will examine this assumption below). 
In addition, the spin flip-flop term includes broadening 
of the inhomogeneous linewidth Tg of the spin transition 
with magnetic field, i.e. Fg = F^-l-yB. This also leads to 
the anti-hole broadening discussed in the previous anal¬ 
ysis, from which we extracted 7 = 6 x 25 GHz/T. 

The second term describes the direct coupling between 
Er ions and resonant, thermally-driven vibrational or 
TLS modes of the glass matrix. While this effect is in 
some ways analogous to the coupling with phonons in 
crystals, one cannot expect it to have the same depen¬ 
dence on the magnetic field and temperature due to a dif¬ 
ferent density of states of vibrational modes in the amor¬ 
phous medium |28| . We therefore include free parameters 
I and m in Eq. 

The third term describes the process where the Er ion 
relaxes through higher-energy vibrational or TLS modes 
of the glass via a second-order Raman-type interaction. 
It is equivalent to the phonon mediated inelastic Raman 
relaxation in crystals, which involves two off-resonant vi¬ 
brational modes. In the case of a crystal doped with a 
Kramers ion such as erbium, the Raman process scales as 
T® [53]; however, since the density of states of vibrational 
modes might be very different in the fiber, the temper¬ 
ature dependence of the Raman process is described by 
the free parameter n. 

By iteratively fitting Eq. to several subsets of ex¬ 
perimental data, we obtain the single set of parameters 
given in Table jlj that allows describing all experimental 
data. The good agreement is exemplified with the two 
datasets shown in Figs. and We note that the value 



FIG. 3: Spin relaxation rate 1/Ta as a function of temperatnre 
at B — 500 G. The dashed lines correspond to the different 
terms of Eq. [^ the Er-Er coupling (blue dashed), the di¬ 
rect process (green dotted), and the Raman process (orange 
dashed dotted). 


for the inhomogeneous line width at zero field, F^, ex¬ 
ceeds the value for REI-doped crystals by several orders 
of magnitude |5DH55] . While we have no explanation yet 
for why 1=1, we note that it is not possible to attribute 
this dependence to a direct-phonon-type process as ob¬ 
served in crystals because such a process must scale at 
least quadratically [53] with the magnetic field. This 
argument is based on the assumption that the phonon 
density of states is not decreasing with frequency, which 
is correct for silica glass at our temperatures |33|. Fur¬ 
thermore, we find that the fitted value for m, describing 
the temperature-dependent part of the coupling to TLS, 
is consistent with values that allow modeling coherence 
properties in REI-doped glasses niiiiiM], and that a 
difference of n from the value in crystals (i.e. ncrystai=Q) 
has also been observed in nanometer size particles [35j . 

Let us now have a look how the interplay of the three 
different processes in Eq. - represented as individual 
lines in Figs. and [^ - determines the lifetime of the 
spectral hole. Fig. shows that the decay rate increases 
monotonically with temperature, and that, for higher 
temperature, it is dominated by the Raman-type inter¬ 
action. Furthermore, Fig. [^ shows that spin flip-flops. 
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FIG. 4: Dependence of the rate XjTa on the magnetic field at 
T = 0.8 K. The solid line shows the theoretical prediction of 
Eq.[^with the set of parameters in Tablejl] and the dashed and 
dotted lines correspond to the individual terms (first term: 
blue dashed, second term: green dotted, third term: orange 
dashed dotted). 



FIG. 5: (a) Dependence of the rate l/Tf, on the magnetic 
field for fibers with different erbium doping concentrations 
and lengths, (b) Absorption profiles of the three fibers, (c) 
Hole depth as a function of magnetic held. All measurements 
are taken at 0.8 K, at 1532 nm for hbers 1 and 2, and at 
1536 nm for hber 3. Red: hber 1 (190 ppm, 20 m), black: 
hber 2 (200 ppm, 10 m), blue: hber 3 (1200 ppm, 3 m). 


i.e. the exchange of spin states between neighboring, res¬ 
onant erbium ions through magnetic dipole-dipole inter¬ 
actions, dominate the decay rate at small magnetic fields 
and 0.8 K. However, the contribution of this process de¬ 
creases rapidly as the field is increased due to the addi¬ 
tional inhomogeneous broadening described by yH, i.e. 
a decreasing probability for neighboring ions to be reso¬ 
nant. For higher fields, the decay rate is dominated by 
coupling to thermally-driven TLS and vibrational modes. 
Hence, for long-lived holes, it is always beneficial to work 
at the lowest possible temperature, but there is an opti¬ 
mum magnetic field. For our fiber, the best value for 
is obtained at 0.65 K and around 500 G. 

Next, we characterize and compare the magnetic field 
dependence of the long decay rates Tt at T = 0.8 K for 
three fibers of different Er-doping concentration: fiber 
1 (the same as in the investigations described above) - 
190 ppm (INO S/N 404-28252), fiber 2 - 200 ppm (INO 
S/N 402-28254) and fiber 3 - 1200 ppm (INO S/N 502- 
28255). Fig. (b) shows the absorption profiles for the 
three fibers. The profiles for the two low-concentration 
fibers (fiber 1 and 2 ) differ greatly, which we attribute to 
a difference in the co-dopants, and those of fiber 2 and 3, 
which features a very different erbium concentration, are 
comparable, suggesting similar co-dopants. The results 
of the rate measurements are shown in Fig. (a). We 
emphasize that the rates for fibers 1 and 2 , which feature 
almost identical Er ion concentrations, are very similar 
over the entire range of magnetic field strength. Contrary 
to what one might expect, this seems to suggests that the 
spectral hole burning properties are not much affected 
by the co-dopants. However, comparing with the data 
obtained for fiber 3, we find that the decay rates strongly 
depend on the erbium doping concentration. To fit the 
data, we therefore include the concentration dependence 
only into ai,a 2 and as, and we assume that the spin 


inhomogeneous broadening, characterized by F^ and 7 , 
is the same for the three fibers. 

As shown in Table |lj we find that the value for ai in¬ 
creases with Er doping concentration, which is consistent 
with the assumption that the first term in Eq.[ 2 describes 
spin flip flops. The concentration dependence of a 2 sug¬ 
gests that the dominant relaxation mechanism at high 
magnetic fields is through spin-elastic TLS involving mo¬ 
tion of Er ions rather than simple vibrational modes of 
glass - or TLS modes created by Er ions [52] . For conclu¬ 
sive information about the scaling of these parameters, 
more investigations using fibers with other doping con¬ 
centrations are needed. Furthermore, we find that the 
third term does not contribute to the decay rates mea¬ 
sured at 0.8 K. 

Finally, we assess the quality of the spectral hole burn¬ 
ing in terms of the hole depth for these three fibers after 
50 ms waiting time and at T = 0.8 K - the results are 
shown in Fig. (c). We find that the erbium ion con¬ 
centration has a direct impact on the hole depth - the 
smaller the concentration, the deeper the hole, which is 
due to reduced spin flip flops. Furthermore, we see that 
the optimum magnetic field also depends on the doping 
concentration, which is a consequence of the change in 
relative importance of the spin flip-flop process versus 
coupling with TLS. 

In conclusion, we established the existence of 
long-lived, narrow and deep persistent spectral holes 
generated via spin-level storage in a rare-earth-ion 
doped amorphous host material. More precisely, using a 
weakly doped silica fiber under an optimized magnetic 
field, we observed population storage in electronic 
Zeeman levels of erbium ions with lifetimes approaching 
a minute. We furthermore developed a model that 
identifies lifetime-limiting mechanisms, in particular 
spin flip-flops at low magnetic fields and coupling to 
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TLS and vibrational modes at high fields. The model 
shows that the use of an amorphous host, which results 
in large spin inhomogeneous broadening, and low doping 
concentration leads to small spin flip-flop rates even 
at low magnetic fields, and hence to the possibility of 
observing deep and long-lived spectral holes. Our study 
shines light on the fundamental interaction between 
impurities and vibrational modes in glasses, which, 
as opposed to crystals, is still not well understood. 
Furthermore, our findings allow parameter optimiza¬ 
tion in future applications of persistent spectral hole 
burning, including optical quantum memories, photonic 
processors, configurable filters, and long-lived optical 
storage elements for fiber-optic communication. We 
emphasize that our investigations have been carried 
out on the 1532 nm transition in erbium, which opens 
these applications to the convenient telecommunication 
c-band. 
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